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Abstract: We study matrix models in the /3-ensemble by building on the refined recursion 
relation proposed by Chekhov and Eynard. We present explicit results for the first /3-deformed 
corrections in the one-cut and the two-cut cases, as well as two applications to super symmetric 
gauge theories: the calculation of superpotentials in M = 1 gauge theories, and the calculation 
of vevs of surface operators in superconformal M = 2 theories and their Liouville duals. Finally, 
we study the /3-deformation of the Chern-Simons matrix model. Our results indicate that this 
model does not provide an appropriate description of the Q-deformed topological string on the 
resolved conifold, and therefore that the /3-deformation might provide a different generalization 
of topological string theory in toric Calabi-Yau backgrounds. 
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1. Introduction 

Matrix models in the 1/N expansion have become a powerful tool in the study of supersymmetric 
gauge theories and string theories. For example, as shown by Dijkgraaf and Vafa in [19], the 
all-genus free energies of type B topological string theories on certain non-compact Calabi-Yau 
manifolds can be computed from the 1/N expansion of simple, polynomial matrix models, and 
this leads to exact results for the superpotentials of a large class of N = 1 supersymmetric 
theories [20] . Other applications include the matrix model formulation of Chern-Simons theories 
[45, 46] and the matrix model-inspired remodeling of the B-model [47, 8] for the mirrors of general 
toric geometries. As a consequence of these relationships, the recent progress in solving the 1/N 
expansion of matrix models [24, 27] has found many applications in string theory and gauge 
theory. 

Most of these applications involve the standard Hermitian matrix model ensemble. There 
is a well-known one-parameter deformation of this ensemble, usually called the /3-ensemble or 
the /3-deformation, which involves an extra parameter (3. The standard Hermitian ensemble is 
obtained when = 1, and the special values (3 = 2 and j5 = 1/2 correspond to Sp(iV) real 
quaternionic and SO(iV) real symmetric matrices, respectively. The 1/N expansion for the more 
general, /3-deformed ensemble, has been worked out in an algebro-geometric language by Chekhov 
and Eynard [15, 14]. As in [24, 27], explicit expressions for the expansion of correlators and free 
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energies are obtained through a "refined" recursion relation based on the spectral curve of the 
matrix model with /3 = l 1 . 

The general (3 ensemble also has many applications. For example, the special values /3 = 
2,1/2 lead to the enumeration of non-orientable surfaces (see for example [11, 51]), and this 
can be used to construct non-critical unoriented strings in an appropriate double-scaling limit 
[11, 32] (see [18, 52] for a review of these ideas). These ensembles also appear naturally when 
one applies the techniques pioneered by Dijkgraaf and Vafa to supersymmetric gauge theories 
with SO(iV) and Sp(iV) gauge symmetry [43, 42, 3]. More recently, there has been renewed 
interest in the /3-ensemble in the context of the so-called AGT correspondence between N = 2 
gauge theories and Liouville theory [5]. In this correspondence, conformal blocks in Liouville 
theory are identified with O-deformed partition functions [53] in N = 2 theories, and it has been 
argued in [21] that the general ^-deformation of J\f = 2 super conformal field theories can be 
implemented by a /3-deformed matrix model with a Penner-type potential. 

In this paper we analyze in detail the recursive proposal by Chekhov and Eynard and its 
concrete implementation in various examples. In section 2 we thoroughly study the algebro- 
geometric solution of the loop equations for the /3-deformed eigenvalue model; in doing so, we first 
find a correction to the diagrammatic solution of [15], which was also very recently pointed out by 
Chekhov [14], and discuss various technical issues associated to the /3-deformation with respect 
to the ordinary (3 = 1 case. We moreover present explicit formulae for the very first corrections 
to correlators and free energies in the /3-ensemble for a variety of situations and potentials; in the 
one-cut case and for polynomial potentials, some of these formulae were already derived in [30] 
and used there to analyze the universality properties of the asymptotic enumeration of graphs 
in non-orientable surfaces (see also the recent paper [7] for another derivation of explicit one- 
cut formulae). In section 3 we use these results to study applications to supersymmetric gauge 
theories. The first application is the computation of M = 1 superpotentials, where we recover and 
generalize previous results in [43, 42, 3, 35]. Our second application is to the AGT correspondence, 
where we consider surface operators [6, 41, 22] in a very simple example associated to a sphere 
with three punctures. In this case, we generalize the B-model computation in [41] and show that 
the /3-deformed correlators obtained with the formalism of [15] lead to correlation functions in 
Liouville theory for general background charge. 

One motivation for the present work was to find a matrix model formulation of topological 
string theory in an O-background. This background provides a one-parameter deformation of 
topological string theory (at least on certain toric Calabi-Yau manifolds) which was originally 
obtained via a five-dimensional version of Nekrasov's partition function [53]. The Q-deformed 
topological string was reformulated later on in terms of the refined topological vertex [36]. More 
recently, the holomorphic anomaly equation has been generalized to the ^-background for N = 2 
gauge theories [44] and more generally for the A- model on local Calabi Yau manifolds [34], thus 
providing an important step towards a B-model version of this deformed theory. 

It is natural to try to extend the remodeling of the B-model [8] to this deformation, and the 
refined recursion relation of Chekhov and Eynard is a natural candidate for this, as suggested 
by the arguments of [21] and by our computations in Section 3.2. In order to test this idea we 
analyze, in section 4, the /3-deformed Chern-Simons (CS) matrix model of [45]. When (3 = 1 this 
model is dual to type A topological string theory on the resolved conifold, and its /3-deformation 
is a natural candidate for the ^-deformation of this theory. Our explicit computations, verified 
by perturbative calculations, show that the recursion of [15] works perfectly well for the CS 

1 This recursion has been reformulated in [25, 17] in terms of "quantum algebraic curves," but the original 
formulation in [15] is more useful for the purposes of this paper. 
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matrix model 2 , but unfortunately they do not seem compatible with the Sl-deformation, at least 
when taken at face value (this was mentioned as well in [34]). An interesting feature we discover 
is a highly involved analytic dependence of /3-deformed amplitudes on the closed string moduli 
with respect to their "refined" counterpart. This degree of sophistication only increases when 
moving to multi-cut models, where the exact formulae we find, e.g. for the cubic matrix model, 
display a more intricated analytic structure as compared to oriented, open amplitudes at higher 
genus [27, 8]. In particular, they cannot be immediately related to the same type of holomorphic 
quasi- modular forms of the ordinary topological string in a self-dual background [2] , and it would 
be interesting to see what kind of generalization would be needed to encompass this more general 
case. 

Our work indicates that the matrix model /3-deformation can be defined and computed for 
the mirrors of other toric Calabi-Yau manifolds. An important example are the mirrors of A p 
fibrations over P 1 . These models can be described by Chern-Simons matrix models on lens spaces 
[1], and one can generalize the computation performed in section 4 to this more general setting. 
In fact, it is likely, in view of the progress in formulating the /3-deformation in a geometric 
language [17], that the /3-deformation provides a generalization of the B-model for the mirrors 
of toric Calabi-Yaus. According to our explicit results, it seems that this deformation will be 
in general different from the O-deformation. If this is indeed the case it would be interesting to 
understand more aspects of this deformation. For example, one could use the Chekhov-Eynard 
recursion, together with the strategy of [26], to formulate a holomorphic anomaly equation for 
the /3-deformed free energies. More generally, one should try to understand this deformation in 
the language of the A-model and in the gauge theory language. 

2. Beta ensemble and topological recursion 

In this section we review and analyze the formalism of Chekhov and Eynard [15], which pro- 
poses a topological recursion for the beta ensemble of random matrices. We will discuss the 
implementation of their formulae and present explicit expressions for various models. 

2.1 General aspects 

In terms of eigenvalues, the beta ensemble of random matrices is defined by the partition function 



In what follows we will mainly follow the normalizations in [15]. The connected correlators are 
defined through 



The correlator W(p) for h = 1 is usually called the resolvent of the matrix model. Both the free 
energies and the connected correlators have an asymptotic expansion in g s , in which the 't Hooft 
parameters are kept fixed. In the case of the free energy F = log Z, we have 




(2.1) 




h>\. 



(2.2) 




(2.3) 



k,l>0 



2 This is not entirely guaranteed a priori, as the formalism of [15, 14] applies in principle to polynomial or at 
most logarithmic potentials. 
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where 

7 = " a/F 1 - (2.4) 
For the first few terms we find, explicitly, 

F = g; 2 PF 0fi + g-\p - l)F 0i i + (/3 + /T 1 - 2)F , 2 + F 1>0 

+ ^(^^^3 + (1-^1,1) (2 . 5) 

+ 9 2 s (V ^2,0 + {1 ^F 1 , 2 + ^3^^0,4) + • • • 

The g s expansion of the connected correlators is written as 

00 

W( Pl , ..., Ph ) = J2 h 29 W g ( Pl , . . . , Ph ), (2.6) 
9=0 

where g can be an integer or a half-integer, and h is defined as 

h = W (2 - 7 > 

This expansion defines the "genus" g correlators, which can be in turn expanded as 

[9] 

W g { Pl ,..., Ph ) = Y,l 2g ~ 2k W k ,2 9 -2k{pir-- ,Ph), (2-8) 

k=0 

and leads to the following expansion for connected correlators, 

/ \ h ~ ^ 

(TiM ni ■ ■ ■ TvM nh )^ = ( ^ ) Yl Yl h 2 9- 2 7 2 9~ 2k (TrM ni ■ ■ ■ TrM nh )[ c J. (2.9) 

9>0fc=0 

The beta ensemble might be regarded as a natural deformation of the standard Hermitian 
ensemble, since when (3 = 1 (2.1) becomes the standard partition function of the (gauged) 
Hermitian matrix model. In this case, in the expansion of the free energy and the correlators 
only the terms F g> o and W g o contribute (with g a non- negative integer). This leads to the 
standard expansion in powers of g 2 of the Hermitian matrix model. On the other hand, there 
are two special values of /3 which have a matrix model realization: for (3 = 1/2, (2.1) describes 
an ensemble of real symmetric matrices with orthogonal SO(N) symmetry, while the case (3 = 2 
describes an ensemble of quaternionic real matrices with symplectic Sp(iV) symmetry. 

Example 2.1. The Gaussian (3 ensemble. In the Gaussian case 

V{x) = x 2 (2.10) 
the matrix integral (2.1) can be computed at finite N by using Mehta's formula 

/ ndA l |A(A)|^ e -i^ 1 ^ = (2vr)^n • (2-H) 



i=i j 
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The result can be expressed in terms of the double Gamma Barnes function 

r 2 (x|o,6) = exp^| i=o C2(a;a,6,a:)) , (2-12) 
where the r.h.s. involves the Barnes double zeta function 

1 roo -tx 

Ua , aM = W) j a dtt .-i_____, (2 . 13) 

see [56] for a summary of properties of these functions. Indeed, it is easy to show that 

N 

\[T{l + pj) = {2ir) N l 2 p N l 2+ P Ni - N - l V 2 T{l + N p)T(N)T 2 l {N;l/ (2.14) 
i=i 

We can now obtain the large N expansion of (2.1) by using the asymptotic expansion of the 
Barnes double-Gamma function [56], 



, tw ,n 1 / 1 2l 3 2 \ 1 fa b\ , 1 / a b\ , 

logr 2 (x;a,6) = - I --x log* + -x + - ( - + - xlogx - - f 2 + - + - ) log* 



oo 



(2.15) 



- X '(0; a, b) + J^(n - 3)!e„_ 2 (a, b)x 2 ' n , 



n=3 



where e n (a, 6) are defined by the expansion 

1 



(l_e-«t)(i_ e -w, ; _ 2 
and x( s ; &) is the Riemann-Barnes double zeta function, 



= ]T e n (a,b)x n (2.16) 



X(s;a,6)= (am + 6n)" s , (2.17) 

(m,n)eNg 



with Nq = N 2 \{(0, 0)}. Up to some additive terms, one finds 



F = \(it 2 (log(t) - |) g~ 2 , (log(/3t) - 1) g; 1 

+ 12/? l0g(/?t) + + 720/^2 & 



(2.18) 



where as usual 

t = g s N (2.19) 

is the 't Hooft coupling. From this expression we can read off the different of the Gaussian 
ensemble. The asymptotic expansion (2.18) can be written as 

F = -logT 2 (t;-g s ,g s /(3). (2.20) 
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2.2 The Chekhov Eynard recursion for the beta ensemble 

When (3 = 1, the full 1/N expansion (2.6), (2.3) of the matrix model was obtained in [24, 16] in 
terms of residue calculus on the spectral curve of the model. We recall that the spectral curve is 
defined by the following relation 

W (x) = \ (V'(x) - y(x)) , (2.21) 

where Wq(x) is the planar resolvent. In this paper we will be interested in the case of hyperelliptic 
spectral curves. y(x) can be written as 



y{x) = M{x)^/a{x), (2.22) 

where 

2s 

a(x) = H(x-x i ), (2.23) 
i=i 

and thus realizes the plane complex curve F = {(x,y(x)),x € C} as a 2-sheeted cover of the 
complex plane, branched at x = Xf, if p G T, we will denote by p the conjugate point under the 
projection map to the eigenvalue plane 

(x(p),y(p)) = (x(p),-y(p)). (2.24) 

In the following, we will often denote the eigenvalue location as p, therefore writing x(p) = p for 
the uniformization variable. The function M(p) in (2.22) is also called the moment function. In 
matrix models with polynomial potentials M(p) is also a polynomial. If the potential contains 
simple logarithms, as in the Penner model that we will analyze later on, M(p) is rather a rational 
function. In many situations related to topological string theory, M(p) can be written in terms 
of an inverse hyperbolic function [47]. For future use, we will denote by C a contour encircling 
the branch points and the branch cuts between them. 

The Chekhov-Eynard recursion relation, proposed in [15], gives a solution to the 1/N expan- 
sion (2.3), (2.6) in the general /3 ensemble, in terms of period integrals defined on the spectral 
curve (2.22). As we will show in a moment, one important difference between the recursion 
proposed in [24, 27] and the one obtained in [15] is that, in the first case, the recursion can be 
formulated in terms of residues in the branch points of the curve. However, in the recursion 
[15], the expressions for W g with g half-integer involve contour integrals where the integrand has 
branch cuts, and they can not be reduced to residues at the branch points. 

The starting point to derive the recursion relations are the loop equations of the (5 ensemble. 
In the following we will assume that V{p) is a polynomial of degree d. The loop equations have 
been written down explicitly in [25], and they read, with the notations above, 

V'( Pl )W( Pl , ■ ■ ■ , Pk )-U( Pl ,--- ,p k ) 

d 

= 2W(p 1 )W(p 1 , ■■■ ,p k ) + h 2 W(pi,pi, ■ ■ ■ ,p k ) + fry— W(pi,p2, • • • ,p k ) 



dp 



fc-2 



+^1^'^)%-^) (2 - 25) 

3=1 I&K 3 



+ E 



k 

d W(p 2 , • • • ,pj,--- ,Pk)~ W(p 2 , • • • ,Pi, • • • ,Pk) 



dpj Pj - Pi 
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for k > 2, while for k = 1 we have simply 

V'(p)W{p) - U(p) = W 2 (p) + h 2 W(p,p) + hj—W(p). (2.26) 

op 

In these equations, U(p±, • • • ,pk) is a polynomial in p\ of degree 5^2 + deg(U) — 2. It turns out 
that these equations can be solved recursively in the g s expansion. To see this, let us look at the 
simple example of k = 1, and let us plug in the expansion (2.6). The first /3-ensemble correction 
is W 1 / 2 (p). It satisfies the equation 

(V'(p) - 2W (p)) W 1/2 (p) - U 1/2 (p) = 7^^- (2-27) 



This can be solved as, 



7 dW (p) , U 1/2 {p) 



Notice that the r.h.s. in this equation is not a rational function, as it happens in the solution 
of the loop equations in the f3 = 1 case, since the derivative of the planar resolvent involves the 
multivalued function y/a(p). However, one can still use the techniques developed in [24, 16, 28] 
in order to give an explicit expression for Wi/ 2 (p). Let dS(p, q) denote the unique third kind 
differential on the spectral curve having a simple pole at p = q and p = q with residues +1 and 
— 1 respectively and vanishing ^-periods. We can write 

W 1/2 (pi) = -Re Sp=Pl dS( Pl ,p)W l/2 (p), (2.29) 

where p\ is a point outside C. We now take into account that Wi/ 2 (p) has no residues at points 
away from the contour C, as well as no residue at p = oo. The first fact follows from the 
assumption that there are no eigenvalues of the matrix model away from the cut (see [24], eq. 
(2.13)), and the second fact follows from the expansion at infinity expressing Wi/ 2 (p) in terms 
of correlation functions, 



p" 

n>l 1 



By contour deformation, we find that 



Wy 2 (jp) = ^ jf d5 (P> <l)Wi/2(q). (2.31) 
Using now the loop equation and the expression for the spectral curve, we find 

Since Ui/ 2 (q) is a polynomial in q, the last integral vanishes, and we obtain 



(2.33) 



The same result can be obtained using the inversion operator of [16, 15]. Using the expansion 
(2.8), we can rewrite this as 

WaM = -^J ciSM m y (2 .34) 
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where we have assumed that V"{q) is analytic inside C. 

Let us now consider the case k = 2. The first non-trivial correction in the (3 ensemble to the 
two-point function is Wi/ 2 (p\,p 2 ). It satisfies the equation 



(V'( P1 ) - 2W (pi)) W 1/2 ( Pl ,p 2 ) - U 1/2 ( P1 , P2 ) 

d d 
= 2W l/2 (pi)W (p 1 ,p 2 ) + j—W {pi,p 2 ) + 



W l/2 (p 2 ) - W l/2 {pi) 

P2 - Pi 



(2.35) 



We can use the same contour deformation argument. There will not be any contribution from 
the polynomial Ui/ 2 (pi,p 2 ) nor from 



d 
dp 2 



Wi/ 2 (p 2 



P2 - Pi 



However, there is a contribution from 



d 
dp 2 



W 1/2 {pi) 

P2 ~Pl 



W 1/2 (p!) 
(Pi ~ P2) 2 



and the final expression is, 
^i/2(Pi,P2) = f r~\ — 



Ic y(p) 



1 1 \ d 

2W 1/2 ( Pl ) [ W ( Pl , P2 ) + 2 {pi _ p2) 2 ) +lQ p -Wo(pi,P2) 



(2.36) 



(2.37) 



(2.38) 



It involves the "corrected" two-point function as in [24] and subsequent works 3 . Notice again 
that the integrand in the above formula is not a rational function, due to the derivative term. 
One can see that the general solution for the "genus" g correlators W g (pi, ■ ■ ■ ,ph) is 



where 



W g (p, Pl , . . . ,p k ) = f (E E Wh{Q,Pj)Wg-h{q,PK/j) 

+ w g -i(q,q,PK) + 7^vw s _i/ 2 (9,pi,-- - ,Pk)j 

1 $k2<) g 



(2.39) 



Wg(pi,...,P k ) = Wg(pi,...,p k ) + 



(2.40) 



2{p 1 -p 2 ) 2 

The free energies can be computed by using the loop inversion operator introduced in [16, 
15, 27]. In the "stable" case, i.e. for (k,l) / (0,0), (0,1), (1,0) and (0,2), they are given by 

where 



&(q) = y(q) (2.42) 

is a primitive of the spectral curve. For the unstable cases, we have specific formulae which can 
be found in [15, 14]. In this paper we will be particularly interested in the first correction to the 
free energy, which is given by 



*b,i = J dq\y{q)\ log \y{q)\. 



(2.43) 



*This correction does not appear in the formulae of [15], see [14] for a careful statement of the recursion. 
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Here, the integration is over the union of the intervals where the density of eigenvalues is non- 
vanishing. To make our notation simpler, we have denoted this support by C again. 

In order to obtain concrete results for the correlators using (2.39), we need explicit formulae 
for the differential dS(p,q). When the spectral curve is of the form (2.22) we can proceed as 
follows [24] . We define the Aj cycle of this curve as the cycle around the cut 

(x 2j -i,x 2j ), j = l,---,s-l. (2.44) 

There exists a unique set of s — 1 polynomials of degree s — 2, denoted by Lj(p), such that the 
differentials 

w = ^±0=&p (2.45) 

satisfy 

Ui = 5ij, i,j = ,s- 1. (2.46) 



/ 

J A. 



The UiS are called normalized holomorphic differentials. The differential dS(p, q) can then be 
written as 

dS(p,q) = ( ZTH -T. C M)Lj{p) J dp (2.47) 

where 

C 3 (q):=^- i -^=—. (2.48) 
2vri J Aj ^Jcr{p) P-Q 

In this formula, it is assumed that q lies outside the contours Aj. One has to be careful when q 
approaches some branch point Xj. When q lies inside the contour Aj, then one has: 

Cr(q) + = — / (2.49) 

which is analytic in q when q approaches x 2 j-i or x 2 j- 

2.3 One-cut examples 

In the one-cut case we simply have 

dS(p,q) ___ 1 



y(q) M(q)yft(p)(p-q)' 



(2.50) 



We will now present some explicit formulae for the very first corrections to the connected corre- 
lators. 

The first correction to the resolvent is given by (2.51), and we find 

^' w = -^iTOFir (2 - 51> 

An explicit, general formula for this correlator was obtained in [30] by using contour deformation. 
Assuming we have a polynomial potential of degree d, we will write the moment function as 

d-2 

M{z) = cY[{z- Zi ) (2.52) 

i=i 
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where c is a constant. We can calculate (2.33) by deforming the contour. This picks a pole at 
q = p, a pole at infinity, and poles at the zeroes of M(z). A simple computation gives 



W ,i(p) 



1 y'(p) 



+ 



d _ 1 + yV°R 

^ p- Zi 



2 y{p) 2^p) 

This can be written in a way which makes manifest the absence of singularities at p = zf. 



(2.53) 



Wo,i(p) 



d-l 1 2p-a-b 



2 V / ^) 4(p-a)(p-6) 2^p) 



E 



P-Zi 



(2.54) 



In the one-cut case it is also possible to write a very explicit formula for i*b,i (or rather for 
its derivative w.r.t. the 't Hooft parameter t). Using that (see for example [18]) 

2 



9 t y(q) 



we find 



* Jc JWl 



(2.55) 
(2.56) 



jogjj/(g)J 

This is easy to calculate in terms of the parameters (2.52) appearing in the moment function, 
and one finds the general one-cut expression, 

1 / a + b 



1 (b — a\ 2 ^ — v 
d t F , 1 = l + -]og(—j-) +logc + ^log 



+ 



(2.57) 



For higher corrections, general formulae become cumbersome (see [7] for an example), but ex- 
pressions for particular potentials are easy to derive. 



Example 2.2. The Gaussian potential. Let us consider the Gaussian potential, 

„2 



V(x) 



x 
~2 



In this case, the moment function M{p) is trivial and we simply obtain 



1 



Wq(p) _ 

yip) 2 Wp 2 -4t P 2 



p 



At 



(2.58) 



(2.59) 



Higher order correlators can be similarly computed in a straightforward fashion from (2.39). We 
find for example 



W 0l 2(p) 

Wo, 3 (p) 
W hl (p) 

w h2 ( P ) 



p 2 + t 



p 



(p 2 -4tf 2 (p 2 -4t) 2 ' 
/ {p 2 + t) (p 3 + 2pi) ' 



\( p 2-4t) 7/2 {p 2 -Atf)' 
if p 2 + 6t p(p 2 + 30t) \ 

2 ^(p2_ 4t )7/2- ( p 2_4i)4 J' 

1 / 21p 4 + 422p 2 t + 200t 2 3 (7p 3 + 48pt) 



(p2-4t) n / 2 



(p 2 - 4t)" 



(2.60) 
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V(p) 


: / 


i F 


o a 


b 



Figure 1: The cubic matrix model in the single-cut phase. 
Example 2.3. The cubic potential. Let us consider a cubic potential 



V{x) = X - + g -x\ 



(2.61) 



with a classical maximum at p = —1/g and a minimum at p = (see Fig. 1). In the stable 
one-cut phase the eigenvalue density is supported on an interval (a, b) around p = 0: the spectral 
curve then takes the form 



l{p) = M{p)yJ(p- a)(p-b), M(p) = g(p - p Q ) 



(2.62) 



where 



Po 



1 a + b 



(2.63) 



The branch points can be expressed as a function of the (only) 't Hooft parameter by imposing 
the correct asymptotics for the planar resolvent [10], and one finds, as a power series in t, 



a = 2Vi-2gt + Ag 2 ^' 2 - 12g 3 t 2 + 36c/ 4 i 5/2 - 128g 5 t 3 + O (t 7 / 2 ) , 
b = -2\Tt - 2gt - Ag 2 t 3/2 - 12g 3 t 2 - 36gH 5/2 - 128g 5 t 3 + O (t 7/2 



It is now straightforward to compute /3-deformed correlators from (2.39). For example, (2.54) 
gives 



Wb,i(p) 



a(2b -3p + po) - 3bp + bp + 4p 2 - 2pp y/(a - p )(b - p ) -2p + 2p 



4ff(p)(p-po) 



2{p-Po)V a (p) 



(2.64) 
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As an instance we have, up to order t 3 and g 7 



(TrM 3 ) = 



-4gt 3 + 



+ (1 - /T 1 ) (9gt 2 + 118<7 3 t 3 + . . . ) 



+ g s p' 1 (gt + 28g 3 t 2 + 664g 5 t 3 + . . . ) 



f g s {\ - /T 1 ) 2 ( - + VTgt + 182cf r + 2228 5 °r + . . . 



+ # 2 (1 _ /T 1 ) (3^ + 198 5 3 i + 6959g 5 t 2 + 202254 5 7 t 3 + ...) + ... (2.65) 
(TrM 2 TrM) (c) = -4 ( 5 i 2 + 10g 3 t 3 + ...)+ g s (l - /T 1 ) Q + 8 5 t + 106c/ 3 t 2 + 1640 5 5 i 3 ^ + . . . 



Example 2.4. T/ie quartic potential. Consider finally a potential of the form 

V(x) = Y+9X 4 . 

The resolvent is given by [10] 

W Q {z) = \{z + 4gz 3 - (1 + 8ga 2 + 4gz 2 ) \/ z 2 - 4a 2 ) , 
where a is a function of g, t 



_L(_ 1 + N /TTw). 



The moment function has two zeros at 

and (2.54) gives [30] 

1 z 



2 _ 1 + 8ga 2 
z o — ' 



4<7 



W ,i(z) 



+ 



+ 



l l-Aa 2 /z 2 z 2 



2z 2 - 4a 2 2z v / l-ia 2 /z 2 V 1 " Aa ' *> lz<2 z ( z * ~ Z D z2 ~ z < 



2 • 



(2.66) 
(2.67) 

(2.68) 
(2.69) 

(2.70) 
(2.71) 



This expression leads to explicit results for the enumeration of quadrangulations of the projective 
plane MP 2 , see [30] for more details. 

2.4 Two-cut examples 

Let us now consider the two-cut case, where we have s = 2. In this elliptic case there is one 
single integral C\(p) (2.48) to compute, and we can obtain very explicit expressions in terms of 
elliptic integrals [8]: 

2 



Clip) 



Li 



k{p - x 3 )(p - x 2 )y/{xi - x 3 )(x 2 - X4~) 
2 

tt(p - x 3 )(p - x 2 )^J{x 1 - x 3 )(x 2 - x 4 ) 

Ky/( x i - x 3 )(x 2 - X~j) 
2K(k) 



(x 2 - x 3 )n(n 4 , k) + (p- x 2 )K(k) 
(x 3 - x 2 )n(m, k) + (p- x 3 )K(k) 



, (2.72) 
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where 



k 2 = 



(xi - x 2 )(x 3 - x 4 ) 



W 4 



(x 2 - xi)(p- x 3 ) 



(xi - x 3 )(x 2 - x 4 )' (x 3 -xi)(p- 
H(n, k) is the elliptic integral of the third kind, 

U(n,k) = 



X2) 



dt 



(x 4 - x 3 )(p - X 2 ) 
(x 4 - x 2 )(p - x 3 ) 



(1 - nt 2 )y/(l — t 2 )(l — k 2 t 2 



(2.73) 



(2.74) 



and K(k) is the standard elliptic integral of the second kind. The leading correction Wo,i to the 
resolvent is given by (2.51). We will split the r.h.s as W { f (p) + wff(p), where 



1 



y'(g)dg 



4vri v / CT(p) Jc (?)(p - ?) ' 



Cl es (q)L iy i(q)dq , <f C 1 (q)L 1 y'(q)dq 



+ 



(2.75) 



(2.76) 



47ri7^) LJ Cl M(g) J C2 M(g) 

When M(p) is a rational function of p, the integrand in (2.75) is a single valued meromorphic 
function outside the cuts and we can compute W^ip) by deforming the contour and picking 
up poles just as we did for the single cut case. On the other hand, as was pointed out in the 
discussion of Section 2.2, this is not the case for the expressions (2.72) for C\{p) and C[ es (p), 
which are only well-defined in the neighbourhood of the cuts [x3,x 4 ] and [x 4 ,X2] respectively. 
A way to treat the integrals appearing in (2.76) is the following: for a fixed polarization of the 
spectral curve, the elliptic modulus k in (2.73) vanishes by definition when we shrink the vl-cycle. 
By expanding the complete elliptic integrals n(n, k) and K(k) appearing in (2.72) around k = 
and integrating term by term, we obtain an expansion of the form 



W$(p) 



4TTly/a(p) 



E 

n=0 



C^ res (q)L iy '(q)dq 



Ci 



M(q) 



+ 



c f (g)W( g )d g 

M(q) 



where we denoted 



1 d n f 



n 



! dk r - 



k n , (2.77) 



(2.78) 



k=0 



At any fixed order in k, by formulae (A. 2) and (A.l), the integrands of (2.77) are algebraic 
functions of q as long as the moment function is rational, and can be computed exactly in terms 
of complete elliptic integrals. 

It should be stressed that, while (2.77) yields only a perturbative expression valid for small 
k, this procedure holds true for a generic, fixed choice of polarization 4 . It therefore provides a 
way to expand the amplitudes around any boundary point in the moduli space where the spectral 
curve develops a nodal singularity. 

Example 2.5. The cubic matrix model. As a first application of our formulae, let us consider 
the case of the cubic matrix model with 



(2.79) 



4 In particular, it continues to hold true when we vary the choice of A and B cycles, thereby changing the very 
definition of dS(p, q) and k; for example, in the context of Seiberg-Witten curves, this would allow us to find 
expansions in any S-duality frame, also at strong coupling. 
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in the two-cut case. The spectral curve reads 

y(p) = M(p)y/^p), M(p) = g, a(p) = y/(p - Xl )(p - x 2 )(p - x 3 )(p - x 4 ). (2.80) 

Following [13, 39] we can parametrize the branch points in terms of a pair of "B-model" variables 

Oi, z 2 ) as 

^Xj = 2Q, x 2 - xi = 2^/zJ, x 4 - x 3 = 2- v /i 2 ~, -x\ - x 2 + x 3 + x 4 = 21 (2.81) 



where 



Q = —, i = 

9 



9 A 



2{ Zl + z 2 ). 



The 't Hooft parameters can be computed explicitly in terms of complete elliptic intej 
as 

(x 4 - X 3 )(x 2 - X 4 )(xi - X^ 2 



h 
t 2 



^\J{^\ - x 3 )(x 2 - x 4 ) 

(x 4 - X 2 )(x 2 - Xl)(x 3 - X 2 ) 2 



where 



k 



j _ (xi - x 2 )(x 3 - x 4 ) 



II (hi, k) 
II (n 2 , k) 



(xi - x 3 )(x 2 - X 4 )' X 3 -Xl 



x 3 - X 4 _ Xl — x 2 

ni = , n 2 = . 

xi - x 3 



This can be inverted as 



zi = -4ti + \<og 2 t\ - 2Ag 2 t x t 2 + . . 
z 2 = 4t 2 - 2Ag 2 ht 2 + I6g 2 t 2 2 + . . . 



Let us turn to compute W^(p) first. We can write it as 

(b) , , _ h(ti,t 2 ) + h{t\,t 2 ) 



with 



We find 



<Ap) 



n=0 



C [ ? Ucg (q)y'(q) 
Y 4mM(q)L 



k n . 



-h = 32^1 (85 + 24 5 3 z 2 + 8l/z 2 2 ) + (20 5 3 + 208 5 5 z 2 + 1269/ z 2 2 ) + 



/ 2 



5^2 . W z\ 

X + ^6- +Z1 V 4 



3g 3 z 2 I3g 5 z 2 



32 



64 



(2.82) 
$rals [33] 

(2.83) 
(2.84) 

(2.85) 

(2.86) 
(2.87) 

(2.88) 



2 \ + z 2 ( 8l9 b z 2 + 1269g^ ( . 



whereas the residue computation for Wq^ (p) yields 



TX AA)f \ 1 1 1 + 2 9P 



4 X : 
1=1 



(2.89) 



(2.90) 
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We can compare our result to explicit perturbative computations for the /3-deformed cubic matrix 
model, along the lines of [40]. As an example, (2.88), (2.90) together yield up to quadratic order 
in t\ and t 2 

(TrM) = - gt\ - gt\ + 4gt 2 h + 30g 3 t 2 2 h - 30g 3 t 2 t 2 - 708 5 5 t|t? + . . . gj 1 



+ (1 - /T 1 ) [g^! + t 2 ) + g 3 (9*1 - mf) + g 5 (-162£?t 2 - 162^*1) +...] + 



which perfectly agrees with the computation from perturbation theory. 



(2.91) 



Interestingly, a closed form expression for Wo,i(p) can be found as a function of the branch 
points. It was shown in [35] that, for matrix models with constant moment function M(p) = g, 
Wo,i is directly related to the planar resolvent as follows 

W ,i(p) = d t W (p) - jd p \na(p). (2.92) 

The first term of the r.h.s. can be evaluated very explicitly upon expressing the derivative w.r.t. 
the total 't Hooft coupling in terms of derivatives with respect to the branch points, following 
[48]. The partial derivatives Aij = satisfy the linear system 

4 



J]%)^=^ 2l (2.93) 
i=i 

4 

^MixjKiAij = 47n5 i)2 , (2.94) 



i=l 



where we denoted 



K, 



[" V^W. (2.95) 

As the Aij are completely determined by (2.93)-(2.94), it is straightforward to perform explicitly 
the derivatives in (2.92) and obtain a compact expression for Wo,i (p) as a function of the branch 
points. We get 



w 0l i(p) 1 



U(h 2 ,k) tt^(xi - x 3 )(x 2 - x 4 ) 

(X2 " X3) -ff(jb) ^W) — +P ~ X3 



It is worthwhile to remark that this expression has a more involved dependence on the branch 
points as compared to oriented, open string amplitudes at higher genus. The ordinary topological 
recursion [24, 27] prescribes the following general form for the /3 = 1 correlators in the two-cut 
case 

3g-3+2h . . „ 3g-3+2h 

W g , (pi,..., Ph )= E [j+ijr) fn({x i },{P j })= E (Mr)) n fn(r,{Pj}) (2.97) 

n=0 ^ ^ '' n=0 

where r is the half-period ratio on the mirror curve, f n are holomorphic, weight — 2n modular 
forms for fixed pj, and E 2 {t) is the second Eisenstein series (see [9, 12] for a detailed discussion). 
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x 2 


V(p) 


X, 


x 4 


/ 

p 







Figure 2: The double- well potential in the symmetric 2-cut phase. 



In particular, only first- and second-kind elliptic integrals are involved for /3 = 1, whereas in 
the /3-deformed case, as (2.96) shows, we have a more sophisticated dependence on closed string 
moduli due to the appearance of elliptic integrals of the third kind at prescribed values for 
the elliptic characteristic. It would be interesting to track the origin of this higher degree of 
complexity for /3-deformed amplitudes. 

Example 2.6. The symmetric double-well. As the simplest instance of a two-cut model with 
non-trivial moment function, consider the double well potential 

V(p) = -£ + g£, (2.98) 

depicted in Fig. 2. The potential has two minima at p = ±1/J~g and a maximum at p = 0. For 
simplicity we consider the case in which we equally distribute the eigenvalues between the two 
minima, i.e. we restrict to the symmetric slice t± = = i/2, t 2 = 0. The moment function in 
this case takes the form 

M(p) = gp. (2.99) 

The branch points can be readily computed as a function of the total 't Hooft coupling t by 
imposing the Z 2 symmetry between the cuts and the leading asymptotics of the resolvent. We 
get 



xi = -a - + ——, x 2 = -\ — , X3 = \ p-, z 4 = */- + ——. (2.100) 

V 9 y/9 V 9 Va V 9 v 7 ^ V 9 V 9 

We now turn to compute Wq l (p)- In this case, the integrals in (2.76) can be computed exactly. 
To see this, let us consider the PSL(2,C) transformation 

ap + P A = f a P 
7P + 5 ' V 7 5 

with 

a = -5 = (x 4 xi - X2X3) C, 

P = (xix 2 x 3 - X1X3X4 + X2X3X4 - xix 2 x 4 )C, (2.102) 
7 = (xi + x 4 - x 2 - x 3 )C- 



P^P=^4, A=( a ^)ePSL(2,C) (2.101) 
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In these equations, ( is given by 



(■■= 



(x 4 - Xl) 

(xi - Xz)\J (xi - x 4 )(xi - x 2 )(x 2 - x 4 )(x 3 - X 4 ) 



(2.103) 



Let us apply this transformation to the second integral on the r.h.s. of (2.76). Then the sum of 
the two integrals becomes the single definite integral 



y Cl 4m[\M(q) M(q)J 



(xi - x 2 - x 3 + x A )y'(q) 



(xi(x4 - q) + (x 3 - x A )q + x 2 (q - x 3 ))M(q) 



In the case of the symmetric double- well, (2.100) implies 

y'(q) _ y'(g) 



M(q) M(q) ' 



xi + x 4 - x 2 - x 3 = 0, 



therefore 

For Wq^\p) we instead find 

(A) 



<?(p) = 0. 



, , , _ v^4t + 3 V - 3 5 V + 5 (4p 2 + it] 

°' l[p) ~ 2 P9 ^r 



2pg 2 a(p) 



As an example, this yields 



<TrM 4 > 0il = 1 - ^ - t 2 - 2^ 3 - 5 5 2 t 4 - 14 5 3 t 5 + 0( 5 4 ). 



(2.104) 

(2.105) 
(2.106) 

(2.107) 
(2.108) 



3. Applications to supersymmetric gauge theories 
3.1 Superpotentials in Af = 1 gauge theories 

In [19, 20] Dijkgraaf and Vafa argued that superpotentials in a large class of M = 1 supersym- 
metric gauge theories can be computed by using matrix models. Let us consider an N = 1 
supersymmetric gauge theory with gauge group G = U(iV), SO(iV) or Sp(iV), where the super- 
field strength is denoted by W a . There is also a chiral superfield $ in a representation R of 
the gauge group G, with a tree level superpotential W tie e(&), which we will assume to be a 
polynomial of degree d: 

d 

W tTee ($) = J2-™ j - (3-1) 
i=i J 

If all roots of W~l Tee (x) = gd Y\i=i( x ~ a i) are distinct, the matter fields are all massive; a classical 
vev for <I>, where iVj of its eigenvalues are equal to cij, spontaneously breaks part of the gauge 
symmetry, and the massive fields can be integrated out to get an effective action for the unbroken 
gauge degrees of freedom at low energy. 

Depending on the gauge group and the representation we will end up with different patterns 
of gauge symmetry breaking (see the useful summary in eq. (2.1) of [35])). We will be particularly 
interested in the examples where G = SO(N), Sp(iV) and R is, respectively, the symmetric and 
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the antisymmetric representation of the group. In this case, we have the simple patterns of gauge 
symmetry breaking 



so(Ao^n s( w)' 



(=1 

k 



(3.2) 



s P (iv)^n s p(^)' 

where 1 < k < d — 1, and there will be, correspondingly, various gluino superfields for the 
unbroken gauge groups, 

Si = -^Tr(w^ a W^), (3.3) 

where is the superfield strength for the i-th gauge group. According to the proposal of 
[19, 20], the effective superpotential for the glueball superfields, as well as the gauge coupling 
matrix for the infrared-free abelian fields, should be computable from an auxiliary matrix model 
with V(x) = Wt r ee(^)- In particular, the glueball superpotential, as a function of the gluino 
superfields, is given by [37, 43, 35] 

i=i olSi 

where e = ±1 for SO/Sp, respectively and Wyy(Si) is the Veneziano-Yankielowicz superpotential 
(see [35] for a detailed expression). In this equation, -Fo,o> Fo,i are the first two free energies in the 
expansion (2.3), obtained in the /3-ensemble for a matrix model with potential V(x) = Wt re e{x), 
in the fc-cut phase, and with 't Hooft parameters S{. In addition, the gauge theory quantity 

can be computed from the generalized Konishi anomaly [3] and expressed in terms of matrix 
model resolvents [42, 35]: 

T(z) = ±N t d -^l-AeW , 1 (z), (3.6) 

where again e = ±1 for SO/Sp. Similarly, contributions to chiral ring observables induced by a 
non-flat gravity background can be computed in terms of non-planar corrections to the resolvent 

[4]- r , _ , , 

The formulae above for the solution of (2.39) in the polynomial matrix model case give then 

explicit results for computing a large class of vevs of chiral observables for a general Wt re e(%)- 

In particular, Wo,i and .Fo,! yield the unoriented contribution to the the effective superpotential 

(3.4) and gauge theory resolvent (3.6) for a general tree-level superpotential 5 . As an example and 

a test of our computations, let us consider the case of classically unbroken gauge symmetry, where 

k = 1. This corresponds to the one-cut case in the computations above. Using the well-known 

one-cut result (see for example [18]) 

9W ,o (z) 1 



dz " y^(i)' 



(3.7) 



5 In the Appendix A of [35], the expressions for the unoriented contributions to the free energy and the resolvent 
in terms of planar, oriented contributions, are only valid when k (the number of cuts) takes its maximum value 
d — 1 for a given potential. 
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as well as (2.53), we find the general formula 



T(z) 



N 



+ 2e 



2e 



d 



^ P~ Zi 



V^lz) ' y(z) yft@) 
This agrees with the explicit computation for the quartic potential in [3] 

3.2 Penner model and AGT correspondence 

A more sophisticated example is given by the double Penner model 

V(x) = oti log x + «2 log(x — 1) . 



(3.8) 



(3.9) 



This model was recently considered by Dijkgraaf and Vafa 6 [21] in the context of the AGT 
correspondence [5] , where it was shown to give a matrix model representation of the chiral three- 
point function in Liouville theory; its 4d counterpart arises [29] as the dimensional reduction of 
the 6d A\ (2,0) theory compactified on a sphere with three punctures, and is a U(l) theory with 
four hypermultiplets. The spectral curve for this case reads 



p(P -I) 



(3.10) 



where «3 = — t — a\ — «2 and t = g s N. 



An extension of the AGT correspondence in presence of defects was considered in [6] , where 
multiple insertions of surface operators on the 4d gauge theory side were mapped to insertions of 
vertex operators corresponding to degenerate states on the Liouville theory side. In [41, 22] both 
were mapped in turn to ^4-type open topological string amplitudes on the toric geometries that 
engineer the relevant gauge theory. In particular the authors of [41] conjectured and checked that 
the Liouville theory four-point function with one degenerate insertion and vanishing background 
charge Q = b + 1/6 



Z nu ll(P,M 



f\V- a2 /n(l)\ 



h I 



(3.11) 



b=i 



should be expressible in terms of oriented topological string amplitudes computed through the 
Eynard-Orantin recursion applied to (3.10) 



Znuilip, %b 



exp 



l^i\p) + ^4°W) + h (4%) + ^j4°Wj>)) + • • • 



with 



A ( h 9 \ P i, ..., Ph ) = j dpi . . . dp h wjf\ Pl , . . . , Ph ). 



, (3.12) 



(3.13) 



On the other hand, it was proposed in [21] that turning on a background charge Q on the CFT 
side should exactly correspond to the /3-deformation of the matrix model, with the dictionary 
been given by 



3 See also [23] for further developments. 
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This was checked by direct computation in [55, 49, 50, 38] at the level of the free energy. It is 
therefore tempting to look at a combination of the two claims above and compute refined open 
string amplitudes via (2.39), corresponding to degenerate insertions in Liouville theory with non- 
vanishing Q. A natural extension of (3.12) in the /3-deformed case is through an expansion of 
the form 



Fnu\\{p, K b) = log Z nuii (p, h, b) 



E ^ 1 - 9 - k/2 ^---Ph k - (3 - 15) 



n=-l 



g,h,k\2g-2+h+k=n 



The /3-deformed topological recursion allows us to test this proposal in detail. On the CFT side 
it is well-known that Ward identities for the normalized four point function (3.11) reduce to a 
hyper geometric differential equation; more precisely we have that 



Z nul i(p,h,b) =p « (1 -p) ft 2 Fi(A 1 ,A 2 ;B 1 ;p), 



where 



«3 — «2 + Oil 

h 



Ao = b 



Cti — C*2 — 03 

h 



+ Q) , 



2b ai 
IT 



(3.16) 



(3.17) 



By Taylor expanding around H = 0, b = i we obtain 

b 2 



logZ null (p) = jU$(p) + 



+ 



+ 



A 

2b 
4b 2 



b -A^(p,p) + (b 2 + l)Af}(p)) 
(1 + b 2 YAfl{p) - b 2 (b 2 + l)\Afl(p,p) - b 2 A%{p) + b^Afl(p,p,p) 



(b 2 + l) s A<°j(p) - b 2 (l + b 2 )A{%) + b^A^(p,p) - (1 + b 2 
b 2 U { 2 %,p) - b 6 ±Afhp,p,p,p) + h\b 2 + l)Af}(p,p,p) 



+ o 



(3.18) 



On the other hand, we can apply the refined recursion to the spectral curve (3.10). In this 
case, the contour integrals also have contributions from the poles of M(x); as an instance, we 
find for the one-crosscap correction to the resolvent 



2 y{p) 2^p)^ 



1 



where 



The residues give the values 



p-z \M(z) 

Z\ =0, Z 2 = 1, Z3 = OO 

a 2 1 



M'(z) r—- \2z-a-b 



2 v^t 



Ol\ 1 
a p 



a p - 1 ' 



-1, 



, (3.19) 
(3.20) 
(3.21) 
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for i = 1, 2, 3, respectively, and we finally obtain 

W (n) — 1M ' ip) 1 2 P~ a ~ b 1 / ail a 2 1 \ , , 

2M(p) A(p-a)(p-b) 2^{p)\ a oP aop-l)' K ' ' 

The integrated refined amplitudes A^ k (p, • • - p) can be similarly computed in a straightforward 
fashion from (2.39); upon taking into account the dictionary (3.14), we find exact agreement 
with the CFT expansion (3.18). 

4. The /3-deformed Chern— Simons matrix model 

4.1 Definition and relation to the Stieltjes Wigert ensemble 

The /3-deformed Chern-Simons (CS) matrix model on S 3 is defined by the partition function 



^cs(^,^,/3)=^/n^-^^^n( 2sinh ^) 2/3 - 

i=l i<j ^ 



When (3 = 1 we recover the standard CS matrix model considered in [45, 46]. This generalization 
of the CS matrix model is the natural counterpart of the /3-ensemble deformation of the standard 
Hermitian matrix model. 

In [58] Tierz pointed out that the standard CS matrix model could be written in the usual, 
Hermitian form, i.e. with a Vandermonde inteaction among eigenvalues, but with a potential 

V(x)= l -{\ogxf. (4.2) 

This potential defines the so-called Stieltjes- Wigert (SW) matrix model. It is very easy to show 
that (4.1) is, up to a simple multiplicative factor, the partition function of the /3-deformed version 
of the SW matrix model. To do that, we perform the change of variables 

Ui = ce Xi (4.3) 

where c is given by 

c = exp(t-g s (l-(3- 1 )), t = g s N. (4.4) 
A simple computation shows that 

Z CS (N,g s ,(3) = e-^ {logc)2 Z sw (N,g s ,(3) (4.5) 

where 

ZsMN, 9s ,P) = ± [ f[pe-&^<**« )2 X\(u,-u^ (4-6) 
" J i=i 11 i<j 

is the partition function of the /3-deformed SW ensemble. In terms of free energies we have 

F cs (iV, g„P) = -^gj 2 +57 1 (/? - l)t 2 - (13 + /T 1 - 2)- + F SW (N, (3,g s ). (4.7) 
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The change of variables (4.3) has to be taken into account when computing correlation functions 
in the CS matrix model from the SW matrix model, and we have the relationship 



(Tr U n ) sw = exp (nt - ng s {\ - /T 1 )) (Tre nX ) cs 



where 



U = diag(ni, • • • , ujv), X = diag(xi, 



,x N ). 



(4.8) 



(4.9) 



It was shown in [46] that, though both the potential (4.2) and its first derivative are non- 
polyomial, the SW model can be solved at large N with standard saddle-point techniques. In 
particular, the resolvent is given by 



W (p) 

and the spectral curve is 



-log 

P 



1 + e~ f p + ^(1 + e"*p) 2 - Ap 



2p 



(4.10) 



y{p) = M{p)\J {p — a)(p — b) = - tanh 1 



P 



y/(l + e-tp) 2 - 4p 
1 + e~ f p 



where 



M(p) 



. tanh 



-i 



V(l + e-^p) 2 - 4p 



py/(p-a)(p-b) " 1 + e-*p 

and the positions of the endpoints are given by 

a(t) = 2e 2t - e* + 2e%\/f* - 1, 
b(t) 



(4.11) 



(4.12) 



2e 2 * - e 1 - 2e^Ve* - 1. 
For t = 0, a(0) = 6(0) = 1, which is indeed the minimum of (4.2). 
4.2 Corrections to the resolvent and to the free energy 

The SW ensemble is, from many points of view, a conventional one-cut matrix model, and its 
correlation functions and free energies obey the standard recursion relations of [27, 15]. We now 
proceed to calculate the first /3-deformed corrections to the resolvent and the free energy by using 
the recursion of [15]. 

Let us first consider the correction fo the 1-point correlator (2.51). As in the standard 
polynomial case, there is no contribution from V"(q) (since both this function and M(q) are 
analytic on the cut). In order to proceed, it will be useful to change variables from q to ( 
through, 

b - a _ a + b , t ^ t . 

<! = — • (4.14) 

This maps the interval [a,b] to [—1, 1]. Explicitly, 



(4.13) 



In terms of £ we have 

tanh -1 



q = 2eVe*(e* - 1)C + e* (2e* - l) . 



-t n \2 



4q 



1 + e 



-t. 



tanh 1 



C+(l-e-*)" 1/2 



(4.15) 



(4.16) 
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and the moment function reads 

M(Q = 1 - tanrr 1 ( y^~} s TTo ) , ( 4 -!7) 

(ci + c^v^ \C + (1 ~ e-T 1/2 J 

where 

Cl = e 5*/2 v ^3T(2e*-l), C2 = 2e 3t (e* - l) . (4.18) 
The integrand of (2.51) involves then, 

~\ dp M^) = Y q ^ {q ~ a){q ~ b)Aq ~ ^(OdC, (4-19) 

with 

Ve*(e* - 1) ( Ve*(e* - 1)C + e* - l) 

M(0 = ; * 7 t= x ■ (4.20) 

(2Ve*(e*-l)C + 2e* - l) tanh" 1 j^^L^ 



We then obtain, 

1 p-e* 1 e* /" M(C) dC 



Wo l (p) = 7: ~ — 77 , 6 / — 

2 pJa(p) 2p Ja(p)Jcp 



2 P\/<Hpj 2 P y/oip) Jcp- 2eVe'(e* - 1)C - e* (2e* - 1) 2vri ' 



(4.21) 



where the integral involving the first term in (4.19) has been calculated through a contour 
deformation and picking residues at q = 0, oo,p, and the contour C encircles the cut [—1, 1] in 
the C variable. 

We have not been able to calculate the second term in (4.21) in closed form. In order to 
obtain explicit results, we have to perform a series expansion in both p and t. To see an explicit 
example of this procedure, we expand around p = oo to obtain 

Wo,i(p) ^-eV-lJ+e^t) (4.22) 

p -2 

where 

5(t) = - / dC-M(C). (4.23) 

7T1 7-1 

Notice that this integral depends on i only through the variable v = e*. It can be computed 
systematically as a power series in u — 1, which can then be re-expanded as a power series in t. 
We obtain, for the first few orders, 

t 7t 2 71t 3 297U 4 17809t 5 4843t 6 51012187* 7 m . 8N 

5(t) = + (t 8 ) . 4.24) 

w 2 24 720 120960 3628800 5913600 435891456000 K ' y ' 

The planar limit of the vev of Tr U is given by 

(Trc/) 7 = e*(e t -l), (4.25) 
and its first correction is given by (4.22), 

(TrU)IJ = Wb,i(p) = e* {S(t) - e* + l) . (4.26) 
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Together with (4.8) we then deduce that 

CEre*>g? = S(t). 



(4.27) 



On the other hand, a direct perturbative computation of the vev (Tre x ) cs in the CS matrix 
model (4.1) gives 



(Tre*) cs = N + g s 



+ 9s 



+ 9s 



— -1(1- p- l )N 2 - lp~ l N + -(1 - p~ l ) 2 N 
_ 6 24 V 1 ' 2V 8 V H ' 

^1 _ H- N 3 (i - /3- 1 ) - lp~ l N 2 + — (i - /rw 2 

24 720 V M ' 48^ 90 V ^ 7 



(4.28) 



l_ (1 _ r if N + l_ r 2 W _ l)N 



+ 0(g, 



in complete agreement with (4.27). 

Using the same type of techniques we can also compute the first correction to the free energy. 
Using (2.56) we find, 



dp log | M(p) | f b dp log | a (p) \ 



aF SW_,+ f°dp ]ogM(p) f 



J a 2VT 

The last integral can be computed exactly: 

' b dp log \a(p)\ log(l — e~* x ' 



(4.29) 



/ 

J a 



+ 2t. 



47r VWM 2 

The first integral can be written, using again the change of variables (4.14), as 



(4.30) 



dp log \M(p)\ 



1 dC 1 



i * vw 1 



loe 



(C1+C20VW 1 



V 



(4.31) 



where c±, C2 are defined in (4.18). As before, this integral can be computed as a power series in 
t around t = 0. Putting everything together, and taking into account (4.7), we find 



F™(t) = -(log(i) + l)t - - + — + — — — + 0(t 7 ). (4.32) 

2 V &K 1 ' 12 1440 45360 14515200 467775 V ' V ' 



We have again verified the very first coefficients in this expansion against a direct perturbative 
calculation in the CS matrix model. 



It is worth pointing out that the corrections appearing in the CS matrix model when (3 ^ 1 
are much more complicated than the "standard" ones. For example, for f3 = 1 all the correlators 
are polynomials in e*, while the integral giving (4.23) is not. 
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4.3 /3-deformation and the $7 background 

One of the interesting aspects of the conventional CS matrix model with f3 = 1 is that its large 
N expansion equals the 1/N expansion of topological string theory on the resolved conifold [31], 
since it equals the partition function of CS theory on the three-sphere. On the other hand, the 
partition function of topological string theory on the resolved conifold admits a refinement given 
by the K-theoretic version of Nekrasov's partition function for a U(l) theory [54]. This partition 
function can be also obtained from the refined topological vertex of [36]. The first correction to 
the refined free energy of the resolved conifold is simply 



*5j = -Li 2 (e-'). (4.33) 



This expression is much simpler than the result (4.32). 

We can also compare the result for (Tre x ) cs with expectations coming from the theory of 
the refined vertex. When (3 = 1, the correlation function (Tre x ) cs can be expressed in terms of 
the open string amplitude 

ZD(t ' 5s) -2sinh(f)- 

for a D-brane in an external leg of the resolved conifold. Here, Q = e - '. The precise relation 
involves a framing factor, 

/rp X 

The "refined" version of the D-brane amplitude is 7 



{Tre")^ =e t Z D (t,g s ). (4.34) 



r? f. 0\ QVQ2- VOl f, ot-X 

z D {t,g s ,p) = ^ — (4.35) 

92 — 1 



with 



q 



= e~^ 9s , t = e~ 9s/ ^. (4.36) 



Expanding in g s and j3 we obtain 



Z D (t,g s ,/3) = {l-e- t )- + -^=(l-P) + -- - , (4.37) 

so it is clear that the relationship (4.34) is no longer true when we consider the /^-deformed CS 
ensemble in the l.h.s., and the refined amplitude (4.35) in the r.h.s. 

Of course, in the comparisons we have made, we assumed that the 't Hooft parameter tcs i n 
the matrix model is equal to the parameter txs appearing in the refined topological string. We 
have not excluded the possibility that both sides are related by a more general relation of the 
form, 

ics = *ts + /(/Mts) (4.38) 

where /(/3,tTs) vanishes for /3 = 1 (since the two parameters agree in that case). But in order 
to reproduce the above results, the unknown function in (4.38) should be rather complicated. 



7 This form of the amplitude is dictated by an implicit choice of gluing along one of the unpreferred legs of the 
refined topological vertex. Other choices of gluing only result in minor differences in this particular case, which 
for the 0(1) term in (4.37) amount to an overall rescaling by a factor of e*; this obviously leaves unchanged the 
discussion about the comparison with the Chern-Simons matrix model computation. 
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Another possibility is that we have to modify as well the Gaussian potential in order to match 
the fi-deformed topological string. This has been suggested in a closely related context in [57]. 
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A. Useful formulae for elliptic integrals 

In this section we collect a few formulae regarding the expansion of elliptic integrals for small 
values of the elliptic modulus, which are relevant for the computations of Section 2.4. 



K(m) 



oo u 
TT \ -> m 



U(n\m) = - ^ 



2 £ (*!)» 



n k kl 



TIT" 



n 



\ 



7T x - m ' - 



E 

fc=0 



Mi 
(k\) 2 
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